Penetration of a vortex dipole across an interface of Bose— Einstein condensates 
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The dynamics of a vortex dipole in a quasi-two dimensional two-component Bose-Einstein con- 
densate are investigated. A vortex dipole is shown to penetrate the interface between the two 
components when the incident velocity is sufficiently large. A vortex dipole can also disappear or 
disintegrate at the interface depending on its velocity and the interaction parameters. 
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I. INTRODUCTION 

Vortex rings, such as smoke rings blown from smoker's 
mouth, propagate in the direction of the axis of the ring 
with a long lifetime. The two-dimensional (2D) analogue 
of a vortex ring is a vortex dipole, which consists of a 
pair of vortices of opposite circulations, propagating side 
by side. A vortex dipole is a very stable structure and 
carries momentum in the direction of propagation. 

In superfluids, a vortex dipole consists of a pair of 
quantized vortices, whose circulations are quantized to 
±h/m with h being Planck's constant and m being 
an atomic mass. Such quantized vortex dipoles have 
been generated and observed in Bose-Einstein conden- 
sates (BECs) of atomic gases [1-4] and of exciton polari- 
tons [5I, |6|] . A quantized vortex dipole propagates with 
velocity h/{md)^ where d is the distance between vor- 
tices. In a uniform superfluid, d is constant in time and a 
quantized vortex dipole propagates with a constant veloc- 
ity. In an inhomogeneous system, on the other hand, each 
vortex of a vortex dipole follows its own trajectory 
or even remains stationary [3, 7]. 

In the experiment reported in Ref. [2] , a vortex dipole 
was generated in an oblate BEC, which propagated 
through the system. When a vortex dipole reached the 
edge of the condensate, the vortex and antivortex disinte- 
grated and they moved along the edge of the condensate 
in opposite directions. This behavior is similar to that of 
a classical vortex dipole moving towards a rigid wall [8] , 
where the vortex and antivortex disintegrate and move 
along the wall in opposite directions. For both a wall 
and the edge of a condensate, a vortex dipole cannot go 
beyond the boundary. In this paper, we investigate the 
dynamics of a vortex dipole moving towards an interface 
between BECs of different components. Few studies have 
been performed on this kind of vortex dynamics at inter- 
faces even in classical fluids. In classical fluids, dynamics 
of vortex rings moving towards a density interface have 
been studied [ol-fTll. 

In the present paper, we will show that a vortex dipole 
can penetrate the interface in a two-component BEC, 
in which quantized vortices in one component are trans- 
ferred to the other component. A vortex dipole also dis- 
appears or disintegrates at the interface. These behaviors 
depend on the incident velocity of the vortex dipole and 



the atomic scattering lengths that determine the inter- 
facial tension and the width of the interface. For some 
parameters, the cores of the vortex dipole after the pen- 
etration are occupied by the other component. When 
only one vortex is occupied, the occupying fraction tun- 
nels and oscillates between the vortex pair. In a trapping 
potential, a rich variety of vortex dynamics can be ob- 
served. 

This paper is organized as follows. Section [IT] gives 
a basic formalism. Section IIIII numerically demonstrates 
the dynamics of vortex dipoles in an ideal system and dis- 
cusses the parameter dependence of the dynamics. Sec- 
tion IIVI investigates the dynamics of vortex dipoles in 
a trapped BEC. Section \V\ provides conclusions to this 
study. 



II. FORMULATION OF THE PROBLEM 

We consider a two-component BEC in mean-field the- 
ory. The dynamics of the macroscopic wave functions 
ipi and ip2 of components 1 and 2 are described by the 
Gross-Pitaevskii (GP) equation. 
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where rrij is an atomic mass, Vj is an external potential. 
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5-wave scattering lengths between atoms in components 
j and j' {j^j' = 1,2). The wave function is normalized 
as J l^jpdr = Nj, where Nj is the number of atoms in 
component j. 

For simplicity, we restrict ourselves to quasi-2D sys- 
tems in the following analysis. The system is assumed 
to be tightly confined in the z direction by a potential 
Vzj{z) and the wave function is reduced to the form 
^j(r) = (j)j{z)%ljj{x^y)^ where (t)j{z) is the normalized 
wave function of the ground state for the potential Vzj {z). 
Multiplying (l)j{z) to Eq. ([1]) and integrating it with re- 
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spect to z, we have the effective 2D GP equation, 



ih 



dt 



2m2 



V^i +^ii|V^i| 



where is the 2D Laplacian, V^j{x^y) 
Vzj{z), and 
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Qjj' = Gjj' J \(l)j\^\(l)j'\^dz. (3) 

We assume that the interaction parameters Qjj' satisfy 
the immiscible condition, 



^11^22 < ^12- 
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We solve the effective 2D GP equation (|2]) numeri- 
cally. The initial state is the ground state obtained by the 
imaginary-time propagation of Eq. (|2]), i.e., i in the left- 
hand sides of Eq. (j2j) is replaced with —1. The imaginary- 
and real-time propagations are obtained using the pseu- 
dospectral method [13]. The computational size is large 
enough that the boundary condition does not affect the 
results. 



III. DYNAMICS OF VORTEX DIPOLES IN AN 
IDEAL SYSTEM 

We first consider an ideal system to study the dynamics 
of vortex dipoles, where the trapping potential is absent 
and the interface between the two components is straight 
along the y axis. For simplicity, we assume mi = m2 and 
911 = 922 = 9' The width of the interface is characterized 
by the parameter 
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For < A <C 1, the density distribution with the bound- 
ary condition lim^^^-oo '^i = lima^^oo ^2 = is given 
by 
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where no is the density far from the interface and ^ = 
ft/ {mgnoY^'^ is the healing length. The width of the in- 
terface w is thus proportional to ^/V^. The time is 
normalized as t = tVs/S,, where = {gno/mY^'^ is the 
sound velocity. 

A vortex dipole is generated by the method proposed 
in Ref. Q (see Fig. 3 in Ref. [14]). When an attractive 
Gaussian potential produced by a red-detuned laser beam 
is displaced in a quasi-2D BEG, a vortex dipole is created 
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FIG. 1: (Color online) Dynamics of normalized density pro- 
files iV^jl^/no and phase profiles argV^j for gn — ^22 and 
A = gi2/gii — 1 = 10~^. The velocity of the incident vortex 
dipole is Viri — 0. Bit's, where Vs is the sound velocity. The 
direction of the vortex dipole propagation is indicated by the 
arrows in the left panels, and the circulations of vortices are 
indicated by the arrows in the right panels. The time is nor- 
malized as t = tvs/^, where ^ is the healing length. The origin 
of time t = is taken arbitrarily. The field of view of each 
panel is 120^ x 30^. 



in front of the potential and "launched" in the direction 
of the displacement [14]. We create a vortex dipole in 
component 1 using a Gaussian potential with magnitude 
Vq > and width w as 



^±1 = Vg 



-Voexp 



(7) 



The position ro{t) is hnearly displaced between the times 
t = and t = and a vortex dipole is launched from 
the potential towards the direction of ro{T) — ro(0). The 
velocity of a vortex dipole can be controlled by the pa- 
rameters Vq, and the function ro{t). The Gaussian 
potential is located far from the interface, and its motion 
does not affect the dynamics near the interface. The dy- 
namics of a vortex dipole therefore depend only on its 
velocity and not on each parameter in Eq. ([7j). 

Figure [1] demonstrates the typical penetration dynam- 
ics of a vortex dipole through an interface. A vortex 
dipole is generated in component 1 and propagates in 
the -\-x direction, which is perpendicular to the interface. 
We note that the vortex dipole in component 1 is accom- 
panied by a "ghost vortex dipole" in component 2, and 
they are located in the same position (arg i/ji and arg iIj2 
at t = in Fig. [T]). As the vortex dipole propagates 
through the interface region, the ghost vortex dipole in 
component 2 is substantiated (|'^2p/^o at i = 100 in 
Fig. [1]) and then the vortex dipole is completely trans- 
ferred to component 2 (|'02p/^o at i = 200 in Fig. [T]). 
After passing through the interface, the vortex dipole 
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FIG. 2: (Color online) Dynamics of normalized density pro- 
files iV^jl^/no and phase profiles argV^j for A = 0.05 and 
Vin = 0.12i;s. The direction of the vortex dipole propagation 
is indicated by the arrows in the left panels, and the circu- 
lations of vortices are indicated by the arrows in the right 
panels. The field of view of each panel is 120^ x 30^. 



in component 1 becomes a ghost (arg?/^i at ? = 200 in 
Fig.H]). 

Figure [2] shows the case of a thinner interface (larger 
A) and a slower vortex dipole compared with those in 
Fig. [TJ When a vortex dipole approaches the interface, 
the straight interface is curved due to the forward atomic 
flow in front of the vortex dipole (iV^ip/no at i = in 
Fig. [2]). The vortex dipole then significantly deforms the 
interface {{ipjl^/no at i = 140 and i = 160 in Fig. E}, 
which creates a new vortex dipole in component 2. The 
vortex dipole created in component 2 then propagates in 
component 2, whose cores are occupied by component 1 
{i = 300 in Fig. [2]). Generation of such a "coreless vortex 
dipole" by a moving potential is reported in Ref. [12]. 
If the incident velocity of the vortex dipole is slightly 
slower, a greater fraction of component 1 is taken away 
from the interface, which forms an elliptic "bubble" of 
component 1 moving through component 2, as in Fig. 1 
of Ref. 

Figure [3] (a) shows the dynamics for an incident veloc- 
ity Vin slower than that in Fig. [2] with the same A. As the 
vortex dipole approaches the interface, the distance be- 
tween the vortices increases and the interface is upheaved 
[i = 320-360 in Fig. [3] (a)]. When the vortices touch the 
interface, the interface is disturbed and the disturbance 
propagates along the interface [i = 400-440 in Fig. [3] (a)]. 
In these dynamics, vortices are not transferred to com- 
ponent 2, and the vortex penetration as in Figs. [T] and [2] 
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FIG. 3: (Golor online) Dynamics of normalized density pro- 
file iV^il^/no of component 1 for (a) A = 0.05 and Vin = 
0.088vs and (b) A = 0.2 and ^;in = 0.039^;s. The field of view 
of each panel is 60^ x 160^. 

does not occur. 

In Fig.[3](b), the velocity of the vortex dipole is slower 
and the width of the interface is thinner than for Fig. [3] 
(a). Near the interface, the vortex and antivortex sepa- 
rate and move along the interface in opposite directions. 
This behavior of vortices is similar to that of classical 
point vortices near a rigid wall. For an inviscid, incom- 
pressible, and irrotatinal fluid, the trajectories of point 
vortices are given by [8] 

1 1 _ 4 

where d is the distance between incident vortices and 
the wall is located at x = or y = 0. According to 
Eq. (jHj), the distance between the trajectories and the 
wall approaches d/2 asymptotically, which roughly agrees 
with the trajectories of the vortices in Fig. [3] (b). 

A vortex dipole moving toward an interface thus shows 
a variety of dynamics as shown in Figs.[l][3l which depend 
on the incident velocity vin of a vortex dipole and the pa- 
rameter A. Figure m shows the parameter dependence of 
the dynamics of a vortex dipole. The penetration of a 
vortex dipole across the interface occurs for large vin and 
small A. For small Vin and large A, a vortex dipole be- 
haves as if the interface is a rigid boundary. The region 
in which a vortex dipole disappears, disturbing the inter- 
face, is located between these regions. 

The parameter dependence in Fig. [H can be under- 
stood qualitatively. A vortex dipole cannot penetrate the 
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FIG. 4: (Color online) Dependence of the dynamics of a vor- 
tex dipole on A and Vin- The region of "Penetration (empty 
core)" corresponds to the penetration dynamics of a vortex 
dipole as shown in Fig. [T] In the region of "Penetration (occu- 
pied core)" the cores of a vortex dipole after the penetration 
are occupied by component 1, as in Fig. [2] The regions of "In- 
terface disturbance" and "Disintegration" correspond to the 
dynamics shown in Figs. [3] (a) and [3] (b), respectively. The 
circles indicate the parameters used in Figs. [iJSl 



interface, when the energy of the vortex dipole £^vd — 
27rfi^nom~^ log((i/^) [16] is much smaller than the en- 
ergy to deform the interface £^interface ^ o-d, where a 
is the interfacial tension coefficient. Using the expres- 
sion for the interfacial tension coefficient in Refs. |l7l.ll8j. 
a = ^ng^A^/^, which is valid for A ^ 0.3, the inequality 

^vd < ^interface rcduCCS tO 



w <^ d. 



(9) 



where w ~ ^/vA is the width of the interface and ^ d 
is assumed. Equation (j9j) indicates that the penetration 
of a vortex dipole is prohibited for large A and small 
Vyd, which is in agreement with Fig. HI We also expect 
that the penetration of a vortex dipole is prohibited when 
the velocity ^cw of the capillary wave on the interface is 
much faster than the velocity of a vortex dipole ^'vd — 
h/ (md), since the disturbance of the interface spreads out 
rapidly before the vortex dipole is transferred across the 
interface, as shown in Fig. [3] (a). The dispersion relation 
of the capillary wave, cj^ = crA:^/(2mno), gives ^ 
[cr/c/(mno)]^/^, and the inequality Vcw ^ '^vd again leads 
to Eq. (|9|). In fluid mechanics, the ratio of inertial force 
to the surface or interfacial tension force is called the 
Weber number, defined by We = pv'^i/cr^ where p is mass 
density, v is characteristic velocity, and i is characteristic 
length. The substitution of p = mno, v = v^d^ and 
i = d gives We ~ w/d^ and therefore the penetration 
of a vortex dipole is prohibited for We <C 1. 

Figure [5] shows the relation between the incident ve- 
locity and outgoing velocity ^out of a vortex dipole 
before and after the penetration of the interface. For a 
sufficiently large incident velocity, the outgoing velocity 
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FIG. 5: (Color online) Incident velocity Vin of a vortex dipole 
versus outgoing velocity Vout after the penetration across the 
interface for A = 0.01 (circles), 0.05 (squares), and 0.1 (tri- 
angles). The line indicates Vin = Vout- 



exceeds the incident velocity. This is understood from 
the relation between the energy E^^ and velocity '^vd of 
a vortex dipole. 
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i.e., the velocity v^d is larger for smaller energy E^d- 
When passing through the interface, a vortex dipole loses 
energy by disturbing the interface, and therefore the ve- 
locity increases. In Fig. [5l Vout rapidly falls to ^ Vin for 
a small Vi^, which corresponds to the "penetration (oc- 
cupied core)" region in Fig. 21 In this parameter region, 
the energy is given by 
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where mcore is the mass of component 1 occupying the 
cores of a vortex dipole after the penetration. Since mcore 
rapidly increases with a decrease in in this parameter 
region, and since the Vont dependence of the right-hand 
side of Eq. (pT]) is dominated by the second term, Vont 
rapidly falls with a decrease in Vi^- 

We also examined various incident angles and found 
an interesting phenomenon. Figure [6] shows the dynam- 
ics for oblique incidence of a vortex dipole. As shown in 
Figs. [6] (a) and [6] (b), the vortex dipole penetrates the 
interface in a manner similar to Fig. [2l After that, a 
fraction of component 1 occupying the vortex cores of 
component 2 oscillates between the two cores [Fig.[6](c)]. 
This phenomenon is due to the tunneling of a fraction 
of component 1 in an effective double well potential pro- 
duced by the vortex cores of component 2. If the distance 
between a vortex and antivortex of a vortex dipole after 
the penetration is larger, the tunneling rate is smaller and 
no oscillation is observed in the relevant time scale even 
when only one core is occupied by component 1 (data not 
shown) . 



(a) component 1 (b) component 2 
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FIG. 6: (Color online) (a), (b) Snapshots and (c), (d) stro- 
boscopic images of the normalized density profiles IV^^I^ for 
A = 0.05, V — O.lSt's, and the incident angle, tan~^ 0.4. In 

(c) and (d), the images from i = 360 to t = 1200 are su- 
perimposed at time intervals of At = 40. The images in the 
dashed circles in (a) and (b) correspond to those in (c) and 

(d) , respectively. The direction of vortex dipole propagation 
is indicated by the arrows. The field of view of each panel is 
50^ X 120^ in (a) and (b), and 120^ x 80^ in (c) and (d). 



IV. DYNAMICS OF VORTEX DIPOLES IN A 
TRAPPED SYSTEM 

We next study the dynamics of a two-component 
BEC confined in a quasi-2D axisymmetric harmonic trap 
mjUj'^{x^ + 7/^)/2. The system is tightly confined in the 
z direction by a harmonic potential mjUj\z^ {2^ and the 
effective interaction coefficient in Eq. ([3]) has the form 
9n' — [^2/(27ra;)]-^/^Gjj/. We assume that the trap fre- 
quencies are the same for both components and that the 
gravitational sag is negligible. In the following simula- 
tions, we use = 27r X 25 Hz and cj^ = 27r x 1.25 kHz. 
We employ the \F = I^tuf = —1) state of ^''Rb for 
component 1 and \F = 2,mF = —2) state of ^^Rb for 
component 2, where an = 99aB and ai2 = 213aB with 
ae b eing the Bohr radius. In the experiment reported in 
Ref. [191, the scattering length of ^^Rb atoms, a22, was 
controlled using the magnetic Feshbach resonance. We 
assume here that a22 is tuned to a22 = 250aB > an, 
and component 2 surrounds component 1 in the ground 
state. The total number of atoms is 4 x 10^ with an equal 
population in each component. 

We create a vortex dipole by the same method [l^ 
as in the previous section, i.e., the attractive Gaussian 
potential in Eq. ([7]) is applied to the inner component 
(component 1). The Gaussian potential with intensity 
Vq = —200huj and waist w = 1.09 /im is moved as 

^ m j Xo-ut {0<t<T) 

^o(^) = \Xo-^xT {t>n ^ ^ 
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FIG. 7: (Color online) Snapshots of the density profiles \ipj\'^ 
of the dynamics for T — 1.0 ms. The inset magnifies the phase 
profile arg'02 in the dashed square. The arrows in the inset 
indicate the directions of circulation. The unit of density is 
8.4 X 10^^ cm"^ The field of view of each panel is 54.3 x 54.3 
/im. 



with Xq = 4.36 //m, u = 1.1 mm/s, and yo = 0. 
The velocity of a vortex dipole is controlled by varying 
T. We first prepare the ground state for xq = Xq by 
the imaginary-time propagation method [Fig. [7] (a)], and 
then switch to real-time propagation. 

Figure [7] shows the time evolution of the system for 
T = 1.0 ms, for which a vortex dipole with velocity 
c:^ 1.2 mm/s is created and launched in the —x direc- 
tion [Fig. [7] (b)]. The vortex dipole then penetrates the 
interface [Fig. [7] (c)]. After the penetration, the cores 
of the vortex dipole are slightly occupied by component 
1 [see the left panel of Fig. [7] (c)]. When the vortex 
dipole reaches the edge of the condensate [Fig.[7](d)], the 
vortices disintegrate and move in the opposite directions 
along the circular edge (data not shown), as observed in 
the experiment in Ref. [2]. 

Figure [8] shows the case of T = 1.4 ms, which produces 
a vortex dipole with velocity 1.1 mm/s. The vortex 
dipole penetrates the interface [Fig. [8] (a)], and the pair 
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FIG. 8: (Color online) Snapshots of the density profiles 
of the dynamics for T — lA ms. The white arrows indicate 
the directions of propagation. The insets magnify the phase 
profile arg '02 in the dashed squares. The arrows in the insets 
indicate the directions of circulation. The unit of density is 
8.4 X 10^^ cm"^ The field of view of each panel is 54.3 x 54.3 
/im. 



of vortices merges, in which component 1 contained in 
the cores also merges to become a droplet [Fig. [8] (b)]. 
The droplet then turns back to the center, forming a 
vortex-antivortex pair in component 2 [Fig. [8] (c)]. The 
directions of circulation of the vortex pair in Fig.[8](c) are 
opposite to those in Fig. [8] (a). When the droplet reaches 
the inner component, it eventually disappears and the 
interfacial wave remains [Fig. [8] (d)]. 

The behaviors of the vortices in Fig. [8] can be under- 
stood as follows. The fraction of component 1 dragged 
into component 2 experiences a force towards the center. 
When the vortex dipole moves outwards [Fig. [8] (a)], it 
experiences a force in the direction opposite to the prop- 
agation, and the pair of vortices approach each other due 
to the Magnus effect (see Fig. 9 of Ref. pj]). When the 
droplet is pushed [Fig.[8](c)], the flow around the droplet 
forms a vortex dipole. If the droplet propagated further 
without reaching the interface, it would split into two 
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FIG. 9: (Color online) Snapshots of the density profiles 
of the dynamics for T — 2.4 ms. The white arrow indicates 
the direction of propagation. The inset magnifies the phase 
profile arg 7/^2 in the dashed square. The arrows in the insets 
indicate the directions of circulation. The unit of density is 
8.4 X 10^^ cm"2. The field of view of each panel is 54.3 x 54.3 
/im. 

droplets (as in Fig. 2 of Ref. [Tsj). 

Figure [9] shows more complicated dynamics, where a 
vortex dipole with velocity 0.9 mm/s is produced with 
T = 2.4 ms. In this case, the penetration does not occur 
and the vortex dipole disappears at the interface. The in- 
terface is disturbed and the disturbance propagates along 
the interface [Figs. [9] (b) and[9](c)]. Interestingly, the 
disturbance refocuses at the opposite side of the circular 
interface and a vortex dipole is created in component 2, 
which propagates in the +x direction [Fig. [9] (d)]. 

V. CONCLUSIONS 

We have investigated the dynamics of quantized vortex 
dipoles in phase-separated two-component BECs. Solv- 
ing the CP equation numerically, we found that a vortex 
dipole can penetrate an interface between the two compo- 
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nents (say, from component 1 to component 2), in which 
quantized vortices in component 1 are transferred to com- 
ponent 2 (Figs. [T] and [2]). The cores of the transmitted 
vortex dipole in component 2 are almost empty (Fig. [T]) 
or occupied by component 1 (Fig. [2]). When the incident 
velocity is slow or the width of the interface is thin, a vor- 
tex dipole cannot penetrate the interface and the vortex 
dipole disappears, disturbing the interface [Fig. [3] (a)], or 
the vortex and antivortex disintegrate and move along 
the interface [Fig. [3] (b)]. Through systematic numeri- 
cal simulations, we obtained the parameter dependence 
of the dynamics (Fig. 2]) and the relation between inci- 
dent and outgoing velocities (Fig. [5]). When a vortex 
dipole penetrates the interface at an oblique angle, the 
cores of the vortex dipole in component 2 are occupied 
by component 1 asymmetrically, followed by oscillation 
of component 1 between the cores due to tunneling. 

We also found a variety of dynamics for a trapped sys- 
tem (Figs. [7]l9]). After the penetration of the interface. 



vortices of a vortex dipole disintegrate and move around 
(Fig. [7]), or change to a bubble, falling back to the in- 
ner component (Fig. [8]). When a vortex dipole cannot 
penetrate the interface, it disturbs the interface and dis- 
appears. In some cases, the disturbance of the circular in- 
terface focuses at the opposite side and reproduces a vor- 
tex dipole (Fig. [9]). These predicted phenomena can be 
observed in, e.g., a Feshbach-controlled ^^Rb-^'^Rb BEC 
confined in a tight pancake-shaped trap. 
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